Noise induced rupture process: Phase boundary and universal scaling of waiting time 
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A bundle of fibers has been considered here as a model for composite materials, where breaking 
of the fibers occur due to a combined influence of applied load (stress) and external noise. Through 
numerical simulation and a mean-field calculation we show that there exists a robust phase boundary 
between continuous (no waiting time) and intermittent fracturing regimes. In the intermittent 
regime, throughout the entire rupture process avalanches of different sizes are produced and there 
is a waiting time between two consecutive avalanches. The statistics of waiting times follows a 
Gamma distribution and the avalanche distribution shows power law scaling - similar to what have 
been observed in case of earthquake events. This suggests that noise induced intermittent rupture 
process might be the true origin of typical scaling in waiting times and avalanches, independent of 
the length scale involved. 

PACS numbers: 02.50.-r, 05.40.-a, 91.30.-f 



Rupture and breakdown [l|, Q are complex processes 
that occur both in micro and macro scales. Natural 
rupture phenomena like earthquake, land-slide, mine- 
collapse, snow-avalanches often appear catastrophic to 
human society. It is therefore a fundamental challenge 
to understand the underlying rupture process so that the 
losses in terms of properties and lives can be minimised 
by providing early alarms. The same crisis persists in 
construction engineering and material industry where de- 
tail knowledge of the strength of the materials and their 
failure properties, are essential. But the physical pro- 
cesses which initiate rupture, help its growth and finally 
results in breakdown, are not completely understood yet. 

Fiber bundle model (FBM) has become a useful tool 
for studying rupture and failure 0] of composite materi- 
als under different loading conditions. The simple geom- 
etry of the model and clear-cut load-sharing rules allow 
to achieve analytic solutions to an extent that is 

not possible in any of the fracture models studied so far 
by the fracture community. FBM was introduced first 
in connection with textile engineering Q and recently 
physicists took interest in it, mainly to explore the critical 
failure dynamics and avalanche phenomena in this model 
Not only the classical fracture-failure (stress- 
induced) in composites, FBM has been used successfully 
for studying noise-induced (fatigue) failure [ll| - [l5| where 
a fixed load is applied on the system and external noise 
triggers the failure of elements. 

In this Letter we introduce the concept of waiting time 
for a noise induced intermittent fracturing process in 
composite materials under fixed external loading. The 
waiting time is defined as the time (Monte-Carlo steps) 
between two consecutive avalanches in the avalanche time 
series for the entire failure process. Through a mean-field 
calculation we show that in the stress-noise space, there 
exists a robust phase boundary between continuous (no 



waiting time) and intermittent fracturing regimes and 
that can be verified by numerical simulations. In the in- 
termittent fracturing regime we study the distributions 
of avalanches and waiting times for different type of fiber 
strength distributions. Finally to compare our model re- 
sults with real rupture situation at large scale we analyse 
earthquake data series (California catalog). 

We consider first a bundle of N parallel fibers - and a 
load {W = <tN) is applied on the bundle. The fibers have 
different individual strengths (x) which are drawn from 
a probability distribution and the bundle has a critical 
strength a c so that without any noise, the bundle 
does not fail completely for stress a < er c , but it fails 
immediately for a > a c . We now assume that each fiber 
having strength Xi has a finite probability P(a, T) of fail- 
ure at any stress a induced by a nonzero noise T: 
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Here P(<7, T) increases as T increases and for a fixed value 
of T and er c , as we increase a, the bundle breaks more 
rapidly. We simulate this failure phenomenon following 
Eq. (JXJ) in discrete time t. After each failure (at the fixed 
stress a) the total load Na is redistributed among the 
remaining fibers equally and we check at time t + 1, if 
the present stress a(t + 1) = W/N(t + 1) can induce any 
further failure following Eq. ([I]). When the value of a is 
considerably large, it so happens that at every time step 
at least a single fiber breaks until the complete collapse 
of the bundle. This is a single avalanche and there is no 
waiting time fl5| . But as we decrease the initial value of 
a, at a limiting value, in a particular time step t not a 
single fiber breaks. We consider this as a single waiting 
time (tw = 1) and the limiting value of a, at which 
the waiting time appears for the first time is denoted by 
(To- This is the onset of intermittent fracturing process. 
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FIG. 1: Phase boundary (<7n vs. T plot) for three different 
type of fiber strength distributions with N = 20000. Data 
points are simulation results and solid lines are analytic esti- 
mates (Eqs. 3,4) based on mean-field arguments. 



After one waiting time, again another avalanche starts 
and eventually all the fiber break after such finite number 
of avalanches. The number of fibers broken during a 
single avalanche is counted as the avalanche size (to). It 
is obvious that as we increase the value of T, the value 
of cr decreases. When the noise is large, the applied 
load has to be smaller for the emergence of a waiting 
time. Thus stress (a) and noise (T) values determine 
whether the system is in continuous rupture phase or in 
the intermittent rupture phase. 

To determine the phase boundary we can give a mean- 
field argument that at a = <jq, at least one fiber must 
break to trigger the continuous fracturing process. Af- 
ter this single failure the load has to be redistributed on 
the intact fibers and the effective stress must be more 
than do - which in turn enhances failure probability for 
all the intact fibers. Therefore in case of homogeneous 
bundle where all the fibers have identical strength Xi = 1 
(therefore <r c = 1), at the phase boundary NP(ao, T) > 1 
giving 



TV cxp 
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which gives 



Co > 
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In absence of noise T, ao = 1 = a c , which is consistent 
with the static FBM results This analytic estimate 
coincides with the data obtained from simulation (Fig. [I]). 
It shows a nice phase boundary between the continuous 
and intermittent fracturing regimes. 

For heterogeneous cases where fibers have different 
strength and the whole bundle has a critical strength <r c , 



we make the conjecture that 



0o > 



l-Tlog(l/N)' 



(4) 



keeping in mind that in absence of noise T, ctq = a c . 
To verify our conjecture we choose heterogeneous bun- 
dles of N fibers where strength of the fibers are drawn 
from a statistical distribution. We have considered two 
different kinds of fiber strength distributions: (1) uni- 
form distribution of fiber strength having cumulative 
form Q(x) = x for < x < 1 and (2) Weibull distri- 
bution Q(x) = 1 — exp(~x k ) where k is the Weibull in- 
dex (we have taken k = 2.0 and 5.0). Each fiber has 
a finite probability P(cr, T) of failure at any stress a in- 
duced by a nonzero T as mentioned before. Similar to 
the homogeneous case, for a particular value of T, be- 
low a certain value of a, the waiting time appears here. 
One can see that the theoretical estimate of phase bound- 
ary agree with the numerical data for the heterogeneous 
cases (Fig. [IJ. However this agreement was much better 
for homogeneous case. This difference can be explained 
through the amount of randomness involved in the re- 
spective systems. In case of homogeneous bundle there 
is no randomness in the fiber strength - the only random- 
ness is coming from the noise term. Whereas in case of 
heterogeneous bundles - there are two sources of random- 
ness - in the fiber strengths and in the noise term. 

Existence of such a phase boundary has important 
consequences on fracturing study in material failure and 
other fracture-breakdown phenomena. In real situations 
of material/rock fracturing, acoustic emission measure- 
ments can show clearly whether an ongoing fracturing 
process belongs to continuous or intermittent fracturing 
phase. Acoustic emissions [li} are basically sound waves 
produced during micro-crack opening within the mate- 
rial body due to external stress and noise factors. Once a 
system enters into continuous fracturing phase the break- 
down must be imminent. Thus the identification of rup- 
ture phase can predict the fate of a system correctly. 

In the intermittent fracturing phase avalanches of dif- 
ferent size are produced separated by waiting times (tw) 
of different magnitude. This happens for a stress value 
a below ao at a certain noise (T) level. We have studied 
the waiting time distribution for both homogeneous and 
heterogeneous bundles with N — 20000. Each curve can 
be fitted with a Gamma distribution 



D(t w ) oc exp(-t w /a)/t w 
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where 7 = 0.15 for homogeneous case and 7 = 0.26 for 
heterogeneous cases (Fig. [2]). The value of a is the mea- 
sure of the extent of the power law regime and it has 
different value for different type of strength distribution. 
We have also studied the waiting time distribution for 
a fixed value of N, but different sets of values of T and 
a, all of which shows Gamma distribution of the form 
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FIG. 2: The simulation results for the waiting time distri- 
butions for three different type of fiber strength distributions 
with TV = 20000. All the curves can be fitted with the Gamma 
form exp(—tw/a)/tl^ J (dashed line) where 7 = 0.15 for ho- 
mogeneous case and 7 = 0.26 for uniform and Weibull distri- 
butions. In the inset we show the data collapse of the waiting 
time distributions with system sizes for uniform distribution. 



of Eq. [51 For a fixed value of N and T as a decreases, 
the power law region extends longer and thus the value 
of a increases, but the exponent of power law decay re- 
mains same. Again for a certain value of N and a as T 
decreases, the value of a increases without any change 
in the power law exponent. These results imply that the 
power law exponent remains unchanged with variation of 
cr, T and N . 

In the waiting time distributions, power law part domi- 
nates for small tw values and exponential law dominates 
for bigger tw values. The noise-induced rupture pro- 
cess, modeled here, has two basic ingredients, external 
stress a and noise T. The noise term triggers initial rup- 
ture which induces one or more load-redistribution cycles 
that finally enhances the effective stress level on the sys- 
tem. Therefore the initial phase of the rupture process is 
dominated by noise term and as the rupture process goes 
on stress factor becomes more dominating. At the final 
stage the stress redistribution mechanism drives the sys- 
tem toward complete collapse through a big avalanche. 
The inherent global load sharing nature is responsible for 
the power law part of the Gamma distribution -as power 
law usuallycomes from a long range cooperative mecha- 
nism 0,0, The exponential part of the Gamma distri- 
bution is contributed by the noise induced failure factor 
P(cr, T). For large tw values one can eventually treat the 
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FIG. 3: Gamma-fitting (dotted lines) to the waiting time dis- 
tributions in California catalog (1984-2002). 



failures to be independent. If P indicates the noise in- 
duced failure probability within tw then the probability 
D{t w ) = A(l - PfwN _ exp(-Pt w N), where A is a 
constant. The normalization of D(tw) requires A ~ N. 
Though for smaller values of tw, one can not ignore the 
correlations between successive failures (responsible for 
the power law part in D(tw)), the exponential scaling be- 
havior for D(tw) can be easily obtained from the above. 
As shown in the inset of Fig. f2J the plot of D(tw)/N 
against twN gives good data collase for different val- 
ues. Such a data collapse indicates the robustness of the 
Gamma function form. 

Our model for noise induced rupture process is not lim- 
ited to any particular system, rather it is a general ap- 
proach and can model more complex situations like rup- 
ture driven earthquakes. There are evidences of stress re- 
distribution and stress- localisation around fracture/fault 
lines in a active seismic-zone and several factors that can 
help rupture evolution are friction, plasticity, fluid migra- 
tion, spatial heterogeneities, chemical reactions etc. In 
our model such stress redistribution/localisation can be 
taken into account through a proper load sharing scheme 
and noise term (T) can represent the combined effect of 
other factors. To compare the waiting time results of our 
model system with real data, we have anal yse d Califor- 



nia earthquake catalog from 1984 to 2002 [17j. We are 
particularly interested in the statistics of waiting times 
[1814201 ] between earthquake events. First, we set a cutoff 
(m c ) in the earthquake magnitude - so that all earth- 
quake events above this cutoff magnitude will be consid- 
ered for the analysis. The distribution of waiting times 
shows similar variation for different cutoff values. We get 
excellent fitting to data points for all the data sets with 
a Gamma distribution 1181 : 



D(t w ) oc exp(-t w /a)/t w 7) ; 



(6) 
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FIG. 4: Numerical data for avalanche size distributions for 
three different type of fiber threshold distributions with N = 
20000. The straight line has a slope —2.5. 



must be answered to develop a prediction scheme based 
on available precursors prior to failure/breakdown. 

Finally, if the avalanche and waiting time data for 
noise-induced fracturing show some similarities with the 
earthquake time series - then we can analyse the syn- 
thetic time series generated from this model more in- 
tensively - compared to the earthquake catalog - where 
data sets are limited and sometimes not large enough 
to perform a good statistical analysis. The new find- 
ings/tools/concepts from such model analysis can be ap- 
plied in earthquake catalog analysis to explore new areas 
- which has not been been done before. 
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with same 7 (~ 0.1) and different a values for different 
cutoff levels: a = 30, 120, 500, 2000 respectively for m c = 
2.5,3.0,3.5,4.0 (see Fig. [3]). 

In general, avalanches or bursts bear important infor- 
mation of the dynamics of intermittent processes. In 
our model the noise T triggers a rupture process which 
continues through load (or stress) redistribution mecha- 
nism. The avalanche size distributions follow an universal 
power law (D(s) <~ s~^) scaling with exponent £ = 2.5. 
This result (FigQJ demands that such intermittent rup- 
ture process belong to the quasi-static fracturing class, 
where the universality of the exponent value has already 
been established 0. 

Identification of phase boundary is crucial for any dy- 
namical system because a system usually changes its be- 
havior as it moves from one phase to another. As we can 
see in our model, there is no waiting time above the phase 
boundary (continuous rupture phase) and waiting time 
appears below the phase boundary (intermittent phase). 
One can also estimate the failure time of the system ex- 
actly [l5| in the continuous rupture phase. 

The concept of waiting time in such a noise-induced 
rupture process is a new and useful concept which allows 
us to study the avalanche time series with the spirit and 
tools that have been used in earthquake catalog analy- 
sis. The similarities in waiting time statistics and scaling 
forms suggest that slowly driven (noise induced) fractur- 
ing process and earthquake dynamics (stick-slip mecha- 
nism) perhaps have some common origin. In case of frac- 
turing in loaded rocks/materials, such study can help to 
identify reliable precursors which can warn of an immi- 
nent breakdown. We notice, in our model system, magni- 
tude of waiting time reduces gradually towards the break- 
down point. What is the form of this variation? Does 
it depend on the applied stress and noise level? Which 
one is the more sensitive parameter? These questions 
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